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We present an all-optical protocol for detecting population in a selected Rydberg state of alkali 
atoms. The detection scheme is based on the interaction of an ensemble of ultracold atoms with two 
laser pulses: one weak probe pulse which is resonant with the transition between the ground state 
and the first excited state, and a pulse with high intensity which couples the first excited state to 
the selected Rydberg state. We show that by monitoring the absorption signal of the probe laser 
over time, one can deduce the initial population of the Rydberg state. Furthermore, it is shown 
that - for suitable experimental conditions - the dynamical absorption curve contains information 
on the initial coherence between the ground state and the selected Rydberg state. We present the 
results of a proof-of-principle measurement performed on a cold gas of ®^Rb atoms. The method is 
expected to find application in quantum computing protocols based on Rydberg atoms. 

PACS numbers: 32.80.Ee,32.80.Qk,32.80.Rm 


I. INTRODUCTION 

Rydberg atoms coupled to electromagnetic fields form 
a promising system for the physical realization of quan¬ 
tum information protocols []| and quantum simula¬ 
tions 0. In these protocols qubits are realized by a 
set of atomic states, which includes one or potentially 
more Rydberg levels. One requirement of these schemes 
is the ability to measure the Rydberg states’ population 
in order to read out the results of the quantum oper¬ 
ations. For accomplishing this task, most experiments 
with ultracold Rydberg gases use methods including field 
ionization and subsequent detection of electrons and/or 
ions on multi-channel plates or channeltrons Q. These 
techniques offer high sensitivity and - for carefully cho¬ 
sen experimental conditions “ selectivity among the 
Rydberg levels . 

Selective field ionization (SFI) techniques are based on 
the fact that the ionization threshold is different for each 
atomic state, increasing from higher to lower lying lev¬ 
els. Hence, by slowly ramping up the electric field and 
monitoring the electrons/ions over time it is possible to 
deduce the initial populations in each level. However, 
the population of a lower lying Rydberg level cannot be 
probed without destroying the population of any higher 
lying Rydberg state. Therefore this method is not appli¬ 
cable in protocols which require independent probing of 
multiple Rydberg states^. Another inherent property of 
methods based on ionization is that the detected atoms 
are removed from the system and cannot be reused. Al- 
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though this atomic loss is negligible in most cases Q, 
it might be a serious limitation in experiments working 
with only one or a few atoms 0 - 

One alternative to ionization detection methods is all- 
optical probing based on electromagnetically induced 
transparency (FIT) 0. This approach has been suc¬ 
cessfully applied in order to non-destructively probe the 
Rydberg level structure in non-interacting [IJ, [l^ and 
weakly interacting 0 Rydbe rg g ases, as well as in the 
presence of electric fields flSl - fl^ . These experiments, 
however, did not access the population of the Rydberg 
state. On the other hand, an EIT-based scheme for the 
optical detection of Rydberg population HH has been 
proposed and demonstrated in dense atomic clouds where 
the Rydberg blockade allows the spatially resolved detec¬ 
tion of Rydberg atoms. 

Here we propose an all-optical scheme for detecting the 
population in a selected Rydberg state in dilute gases. 
By using a series of laser pulses in FIT configuration this 
technique also allows for distinction between coherent su¬ 
perpositions and statistical mixtures of the ground and 
Rydberg states of the atoms. Since this scheme is based 
on time-resolved observation of the optical response of 
individual atoms it may, in principle, be used down to 
the single atom level. Our method is state selective and 
applicable for testing the population not only in the high¬ 
est Rydberg level of interest (cf. SFI) but any lower lying 
or intermediate Rydberg state. 

We present our theoretical model along with numeri¬ 
cal simulations and demonstrate the scheme in a proof- 
of-principle experiment with a dilute gas of ®^Rb atoms 
showing the detection of the population in an initially 
prepared Rydberg state. Our analysis includes character¬ 
istic effects of Rydberg experiments such as blackbody- 
induced depopulation HIJl, superradiance and 
dipole-dipole interaction |25j|. 
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FIG. 1. (Color online) Illustration of the atom-laser interac¬ 
tion as used for the model of time-resolved electromagneti- 
cally induced transparency (FIT). States |1), |2) and |3) de¬ 
note the ground state, the first excited state and a Rydberg 
state of the atom, respectively. In our experimental setup, 
these states correspond to the 5Si/2(F’ = 2), 5P3/2(F’ = 3) 
and 35 Si/ 2(F' = 2) states of ®^Rb. The atomic transition 
|1) 12) is driven by a weak probe laser with the Rabi fre¬ 

quency flp (red), while the transition |2) o |3) is driven by a 
stronger coupling laser with the Rabi frequency Qc (blue). 7 p 
and 7 c denote coherence decay terms. 5p and 5c are the de¬ 
tunings of each laser to the corresponding atomic resonance. 
The radiative decay from the selected Rydberg state |3) to 
the neighboring states is accounted for by including a reser¬ 
voir state |4). Tij denote the respective incoherent decays, 
consisting of spontaneous emission as well as transitions in¬ 
duced by blackbody radiation. 


II. THEORETICAL MODEL 

Let us consider a cold atomic gas interac ting with two 
laser pulses in an EIT-like configuration [^. One of 
the pulses, resonant with the atomic transition between 
the ground state |1) and the first excited state |2), is 
a weak probe pulse well below the saturation intensity, 
while the other one is a relatively strong coupling pulse 
which is resonant with the atomic transition between the 
first excited state |2) and a selected (arbitrary) Rydberg 
nSi /2 state |3) (see Fig. [T|). A single “reservoir” state |4) 
is used to model the neighboring Rydberg states (l3 |. 

We use a seiniclassical approach for describing the dy¬ 
namics of the system and the laser pulses are taken into 
account through their classical electric field. The atomic 
gas is modeled by a motionless ensemble of atoms. The 
state of the atoms is described by the density matrix 
P = Si j=i Pij I*) 01; where the states |2) and |3) ro¬ 
tate with the atomic transition frequencies a; 2 i and W 23 , 
respectively. 

The time evolution of p{t) is described by the master 
equation 

ihp=[H,p]+U[p\. (1) 

Here the Hamiltonian T-L accounts for the interactions be¬ 
tween the atoms and the laser pulses. The effects due to 


interatomic interactions are considered through dynamic 
effective rates in U[p] along with radiative losses occur¬ 
ring in the system. The Hamiltonian % is written as 

H = -^(Hp|2)(l|+L!c|3) (21-t-h.c.) 

-;i(5p|2)(2|+5e|3)(3|), (2) 

where flp = {Epdi 2 )/h and flc = (^ 0 ^ 23 )/^ are the Rabi 
frequencies of the probe and coupling lasers, with E-p 
and Ec being the electric fields, and di 2 and ^23 the 
dipole matrix elements of the corresponding transitions, 
whereas (5p and Sc are the detunings of the probe and 
coupling laser from the corresponding transitions, respec¬ 
tively (see Fig. [TJ. Although we consider a situation 
where both the coupling and the probe laser are resonant 
with the atomic transitions they drive (i.e. (5p = i5c = 0), 
by including these detunings one can account for poten¬ 
tially uncompensated electric and/or magnetic fields in a 
specific experimental realization. The operator U which 
governs the non-Hamiltonian part of the dynamics reads 
as 


U[p] 


— (2(T13P(T31 

H—^ (2cri2pcr2i 

H—^ (2cr43p(T34 
H—^ (2cri4PCr4i 
+ ^ (2tTiipcrii - 
+ ^ (2o'33pcr33 - 


- 0-33P — P<^33) 

- (T22P — pcr22) 

- CT33P — ^0-33) 

- 0 - 44 /? — ^< 744 ) 
piicr - peril) 
P330’ — po'33) , 


(3) 


where (Tkj = \k) {j\ are the atomic projection operators 
(fc,/e{l,2,3,4}). 

There are multiple sources of non-unitary dynamics in 
the system. One of them is the spontaneous emission 
from the first excited state |2) and the Rydberg state |3) 
which we take into account by introducing radiative de¬ 
cay rates r 2 i and r 32 . Another source, if present, is a 
depopulation of the Rydberg state |3) towards the neigh¬ 
boring Rydberg states. The depopulation may occur due 
to several phenomena depending on the actual realization 
of the system, such as amplified spontaneous emission 
and/or superradiance [2J,[2g as well as induced emission 
and absorption due to the blackbody radiation of the en¬ 
vironment [ 23 , [ 2 ^ ■ Following [ 2 ^ , we take these effects 
into account by modifying the third term in Fq. ([3|) to 
the dynamic effective decay rate 

f34(t) = r34^sp • [p44{t) ■ Psup + 1] + r34_bb, (4) 

with Psup a superradiance parameter, and r 34 _sp and 
r 34 ,bb the effective decay rates caused by spontaneous 
emission and blackbody radiation, respectively. In our 
model we assume that the entire population eventually 
ends up in the ground state |1). This assumption is valid 
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provided the ionization from all involved states is negli¬ 
gible. 

The above mentioned phenomena cause population- 
transfer between the atomic states. In contrast, there are 
a group of processes which do not result in a signihcant 
energy-decay in the system but leads to a relevant coher¬ 
ence loss. One such process is the phase noise of the driv¬ 
ing lasers, which is included into the model through the 
coherence decay rates 7 p and 7 c. Due to the redistribu¬ 
tion of population from |3) to |4), atoms in Rydberg nP^ 
states are present in the cloud at various distances. As 
observed by [1^ , the dipole-dipole interaction with these 
nPj state atoms results in an inhomogeneous broadening 
of the Rydberg nSi /2 state. We take this into account 
by adding an effective dephasing term in Eq. to 7c, 

7c(t) = 7c + 73,dd ■ (5) 

The optical response of the cloud under the effect of the 
two laser pulses is given by the macroscopic polarization 
P = J\fTr[pd] where d = I*) 01 + ^-C.) is the 

atomic dipole operator. The absorption a of the probe 
laser is then given by the imaginary part of the electric 
susceptibility %, 

a{t) = Q(x(t)) = ( 6 ) 

CQfii ip 

where Af is the atom density of the cloud. Note that here 
we make the approximation that the cloud is homoge¬ 
neously irradiated and the propagation effects of the laser 
pulses can be neglected. The time dependent absorption 
signal is thus given by the master equation, which, using 
the operators given in Eqs. @ and ©, reads as 


Pll —2 (^pP 21 — DpPi 2 ) -I- r 2 lP 22 + r 4 iP 44 (7a) 

P 22 =2 (^pPl 2 — DpP 21 — DcP23 + DcP32) 

- r 2 lP 22 + r32P33 (7b) 

P33 =2 (^cP23 — DcP 32) — (r32 + 1 ^ 34)^33 (7c) 

P 2 I =2 [^cP31 ~ ^p (P22 — Pll) + 25pP2l] 

- ^ (r 21 + 7 p)P 21 (7d) 

P31 =2 [^cP21 — DpP32 + 2((5p -|- (5 c)P3i] 

- ^( 7 p + 7c+r32 + f34)p31 (7e) 

P32 =2 [~^pP31 — Dc (p33 — P 22 ) — 2(5cP32] 

P44 =r34P33 — r4ip44 (7g) 

% 1 

Pil =2 {-^pPi2 + 25pP4l) - -(r4i -I- 7 p)p 41 (7h) 



FIG. 2. (Color online) Pulse sequence for the coupling (blue) 
and probe (red) lasers. At t = 0 a fraction of the atomic pop¬ 
ulation is prepared in the Rydberg state. In the simulations 
this fraction is an input value, while in the experiment it is de¬ 
termined by the length of the excitation pulse (Exc.). Before 
this, an optical pumping pulse (Pump) is used to pump the 
atoms to the correct polarization. The time evolution of the 
optical density is monitored after the excitation pulse with 
the probe laser (A, B, C). The coupling laser is added in time 
interval B (EIT pulse). The Rabi frequencies are taken from 
the experimental values in part IIVI 

% 1 

P42 = — 2 (^p7'41 + ^cP43) — 2 ^^21 + r4l) PA2 (71) 
P43 = — 2(^cP42 + 25cP43) 

— 2 (^32 + f34 + r4i -I- 7c^ . (7j) 

When there is no coupling to the Rydberg state (Dc = 0), 
an analytical steady state solution for a can be obtained, 
which will be used for normalization: 

^12 / I -p \ 

Oo —-^ ■ r2i(7p -b r2i) 

eo« 

X ['ylr2i + 2^p{nl + rj,) 

+r2i(452 + 2D2 + r2j]-\ (8) 

III. DETECTION OF THE RYDBERG 
POPULATION: RESULTS OF THE NUMERICAL 
SIMULATION 

We numerically solve the equation system © with the 
pulse sequence of the lasers given by Dp(t) and Dc(t) (see 
the sequence A, B, C in Eig.©. The solution provides a 
description of the population dynamics while the atom is 
being probed by the weak laser on the 11) -(-t- |2) transition 
along with a time-dependent coupling between states |2) 
and |3). Furthermore, through Eq. (© it describes the 
absorption of the probe laser, which we give relative to 
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FIG. 3. (Color online) Dynamics of the populations pkk{t) of 
the atomic states (dotted lines) and the relative absorption 
ftrei(t) (solid lines) of the probe laser induced by the pulse 
sequence A, B, C (see Fig. O for atoms initially prepared in 
the ground state (a), the selected Rydberg state (b) and the 
reservoir state (c). Since the population p 22 is almost zero at 
all times, except for a transient population in the hrst 100 ns 
of B, it is not shown here. 


ao (see Eq. ([5])): 

arei(i) = — (9) 

cto 

In Fig. 13] we show results for the cases where the en¬ 
tire population is initially (a) in the ground state (po = 
p{t = 0) = |1)(1|) or (b) in the selected Rydberg 
state (po = |3) (3|). In (c) the population is split be¬ 
tween the Rydberg state and the reservoir state (po = 
0.7 |3) (3|-b0.3|4) (4|). 

Following the pulse sequence, the time evolution of the 
system can be separated into three major parts. When 
the atoms are initially prepared in the ground state, only 
small changes in the populations are visible (see Fig. 
|3Ka)). However, the reasons for these changes are not 
the same in the different parts of the time evolution. In 
part A, when the atoms are only irradiated by the rela¬ 
tively weak probe laser (Hp <C r 2 i), a small fraction of 
the population is transferred to the first excited state | 2 ) 
by the absorbed light. In part B, the population trans¬ 
fer to state | 2 ) is prevented by the strong coupling laser 
applied on the transition between states |2) and |3) and 


the absorption is reduced, which is the well-known ef¬ 
fect of FIT. The timescale for the transparency to build 
up is defined by the Rabi-frequency flc of the coupling 
laser. If the requirement Hp <C is not fulfilled, the 
transparency is only partial. In this case, the two laser 
fields cause two-photon transitions to state |3), and the 
absorption of the probe pulse is nonzero. This absorption 
level (in case of Jp = i5c = 0 ) depends on Op/Hc and the 
decoherences 7 p and 7 c. A consequence of this effect is 
that different initial populations of Rydberg states cause 
a different absorption level in the equilibrium of part B 
due to the dependence of 7 c on P 44 . In part C, where the 
atom cloud is again only irradiated by the weak probe 
laser, the process is very similar to what happens in part 
A with the exception that there is a small fraction of pop¬ 
ulation is state |3). Consequently, the absorption level of 
the probe laser is smaller, because atoms are missing from 
the ground state. 

In the case of initial population in the selected Rydberg 
state |3), the dynamics of the system only differ in parts 
A and B, provided B is long enough to reach steady state 
FIT. In part A, the absorption of the probe laser is close 
to zero, and slowly increases while a small fraction of the 
population decays from the Rydberg state |3). Since the 
lifetime of the Rydberg states is much longer than 10 ps, 
the amount of population transferred by spontaneous de¬ 
cay is small although not negligible on the ps timescale 
of the pulse sequence. At the beginning of part B the 
population in the selected Rydberg state is transferred 
to the ground state | 1 ) through a resonant transfer from 
state 13) to |2) induced by the coupling laser with Rabi 
frequency He and the consecutive spontaneous emission 
from 12) to |1). For He < r 2 i, this process results in only 
a small increase in P 22 , because a half Rabi cycle induced 
by He between states |2) and |3) would take longer than 
the lifetime of state |2) (see Fig. EKb)). Hence, assum¬ 
ing the initial population is either in the ground state 
|1) or the selected Rydberg state |3), we can determine 
the fractions by monitoring the absorption of the probe 
laser in part A of the time evolution. If there is a way to 
ensure that all the population missing from the ground 
state |1) is in the selected Rydberg state |3), then this is 
indeed sufficient. However, if the probability that a frac¬ 
tion of the population is in another state (for example, 
the interaction scheme to be realized contains more than 
one Rydberg state), the absorption level in part A of the 
time evolution is not enough in itself to give information 
about the population of the selected Rydberg state |3). 

As illustrated in Fig. [3] (b) and (c), the absorption 
of the probe laser in part A is the same for different 
initial states of the atoms as long as the population in the 
ground state is the same. In contrast, the dynamics and 
the equilibrium become significantly different in part B. 
Since only the population in the selected Rydberg state 
is transferred back to the ground state by the coupling 
field, the absorption level in part C also changes with the 
initial population in |3). 

Another result of the simulations is the possibility to 
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FIG. 4. (Color online) Numerically calculated relative ab¬ 
sorption Orel of the probe laser in the beginning of part B for 
atoms prepared in states pm(t = 0) = ^(|1} (1| + |3} (3|) and 
p±{t = 0) = ^(|1) ± |3))((1| ± (3|). The parameters used for 
the calculation are the same as for Fig. [3] with flc < r 2 i in 
(a) and flc ~ r 2 i in (b). This distinct signature of coherent 
states is expected to be experimentally observable and even 
more pronounced if flc is larger. 


obtain information on the initial coherence of the sys¬ 
tem. The time evolution of the relative absorption in 
the beginning of part B for three different initial prepa¬ 
rations of the atoms is shown in Fig. SI These initial 
preparations consist of the same fraction of population 
in the ground state |1) and the selected Rydberg state 
13), but the coherence between these two states is dif¬ 
ferent. One of the initial preparations is the mixed state 
/Om(t = 0) = |(|1) (1| + |3) (31), while the other two prepa¬ 
rations are p±{t = 0) = ■|(|1)±|3))((1|±(3|). Comparing 
numerical calculations for these three cases, we find sig¬ 
nificant changes in the beginning of part B, where the 
FIT did not yet reach equilibrium. The absorption level 
during the rest of the pulse sequence is not sensitive to 
the initial coherence. If the Rabi frequency fic of the cou¬ 
pling laser is on the order of r 23 or higher, oscillations of 
the absorption signal can be observed. 


IV. EXPERIMENTAL SETUP 

For demonstrating the detection of Rydberg popula¬ 
tion with time-resolved EIT, we have conducted an ex¬ 
periment on a cloud of Ri 2 x 10^ ®^Rb atoms at a tem¬ 
perature of « 150 pK. In this experiment the atoms are 
trapped in a magneto-optical trap (MOT), loaded to a 
magnetic quadrupole trap and then released. The time- 



FIG. 5. (Color online) Experimental setup for time-resolved 
EIT measurements. The counter-propagating coupling (blue) 
and probe (red) beams are superimposed on the atom cloud 
and separated with dichroic mirrors. The transmission of the 
probe laser is detected by an avalanche photo diode (APD). 
The intensities of the lasers are controlled with acousto-optic 
modulators (AOM). Both laser frequencies are stabilized us¬ 
ing a frequency comb. 


resolved measurements are started after 3 ms of time of 
flight, in order to ensure that all magnetic fields have 
fully decayed while the effects of atomic motion are still 
negligible. The measurements are performed within 30 ps 
(excitation pulse and probe sequence, cf. Fig. [2|). The 
density and optical density at the center of the cloud, 
measured by absorption imaging, are 5 x 10® cm“^ and 
1.7, respectively. 

The transitions from the ground state |5Si/2,F = 2) 
to the hrst excited state |5P3/2,F = 3) and from there 
to the selected Rydberg state |35Si/2) are driven by 
two lasers with wavelengths of w 780 nm (red, probe) 
and ~ 480nm (blue, coupling), respectively (see Fig. [T]). 
Additionally, we use a repumper to pump atoms from 
|5 Si/ 2,F = 1) back to |5 Si/ 2,A = 2) via |5P3/2,F = 2) 
during the whole pulse sequence. The frequencies of 
both lasers used in the experiment are referenced to a 
frequency comb and controlled with slow servo loops 
(< 100Hz bandwidth). The linewidths of both lasers 
are narrowed to less than 2 tt x 20 kHz with fast locks 
(> 1 MHz bandwidth) to scanning Fabry-Perot interfer¬ 
ometers. As the Fabry-Perot cavities are sensitive to 
acoustic noise, the effective linewidth for the experiment 
can be larger. The red and the blue lasers are aligned in a 
counterpropagating configuration (see Fig. [5]). We use an 
acousto-optic modulator (AOM) in each beam to create 
the intensity envelopes of the pulses. The switching time 
of the AOMs is 50 ns (20 % to 80 % light intensity). The 
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time-dependent intensity of the red laser (probe signal) 
is measured with an avalanche photo detector (Thorlabs 
APD120A/M) and recorded with a digital oscilloscope. 
The time resolution of the setup is 20 ns. 

As shown in Fig. [2l our measurement consists of three 
main parts. First we put all the atomic population 
into the Zeeman-sublevel of the ground state matching 
the polarization of the red laser by optically pumping 
for 14.5 ps. This pulse is also long enough that any 
transient effects due to the switching of the AOM wear 
off before the next pulse. Next, we prepare the ini¬ 
tial state of the atomic cloud with an excitation pulse. 
We always apply the red laser for 2 ps with a Rabi fre¬ 
quency of 27r X 3.3 MHz. If we want no population to be 
transferred to the Rydberg state, the blue laser remains 
switched off and due to the very short lifetime of the first 
excited state, practically the entire population remains in 
the ground state. Applying the blue laser with a Rabi 
frequency of 27r x 1.8 MHz for up to 2ps causes a frac¬ 
tion of the population to be transferred to the Rydberg 
state, with the most atoms being excited in the case of 
the full 2ps pulse. Due to the high Rabi frequency of 
the red laser there is no coherent excitation of Rydberg 
atoms. In the third part (probe sequence in Fig. [2]), we 
use the red laser at a low intensity as the probe laser 
during time intervals A, B and C together with the blue 
laser as the coupling laser during time interval B (FIT 
pulse). 

As a reference, we always add one experimental cycle 
without atoms in order to measure the intensity /ref(0 
of the red laser with the photodetector and another one 
with atoms but no excitation pulse and no coupling pulse 
to normalize the data later on. The parameters are then 
varied from cycle to cycle and the transmitted intensity 
/t (t) of the probe light after passing through the cloud is 
measured. The experiment is repeated several times for 
each set of parameters to reduce photo diode noise. The 
optical density OD is calculated as follows: 


The resulting OD datasets are then normalized by divid¬ 
ing them by the OD dataset that had no excitation pulse 
and no FIT pulse [ODo(<)]: 


ODrel(t) 


OD(t) ^ a{t) 
ODo(t) oq 


( 11 ) 


This relative optical density ODj-ei is comparable to the 
relative absorption a{t)/ao that is calculated in the nu¬ 
merical simulation. To reduce the effect of the acous¬ 
tic noise on our Fabry-Perot cavities, we selected the 30 
datasets where the mean transparency in the FIT pulse 
between 7 ps and 9ps was maximal. We observe that in 
all measurements the relative optical density eventually 
returns to the level before the excitation pulse, ensuring 
that ionization effects are negligible. 

In general this scheme is applicable in situations where 
the optical density can be precisely measured. The reso- 
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FIG. 6. (Color online) Measured optical density for the 
dynamics during the probe sequence. The shaded areas 
are 95 % confidence intervals for the relative optical densi¬ 
ties ODrei(t) obtained from the measurements by applying 
Eqs. (HOI) and Hill for three durations of the excitation pulse 
(Ops for I, 0.5ps for II and 2ps for III). Solid lines repre¬ 
sent £t results for the simulated relative absorption a{t)/ao. 
The reversed order of the absorption signals in part B can be 
explained by dipole-dipole interactions. 



lution is limited by technical noise from the photodiode 
and the digital resolution of the subsequent data acqui¬ 
sition system. For low optical densities it is necessary 
to detect not only the absorption of the atoms, but to 
resolve the FIT signal as well. At high optical densities 
the constant resolution of the intensity measurement ad¬ 
ditionally leads to a lower resolution of the optical density 
due to the logarithmic scaling in Fq. m- For sufficient 
averaging, we estimate the presented scheme to be ap¬ 
plicable in the range of optical densities between « O.I 
and w 4. In principal, the optical density can be lowered 
by detuning the probe laser while maintaining the two 
photon resonance for the FIT condition, which on the 
other hand decreases the contrast for the FIT signal and 
therefore only allows for a limited extension of the range. 


V. DETECTION OF THE INITIAL 
POPULATION: EXPERIMENTAL RESULTS 

We demonstrate the optical detection of population for 
three different initial Rydberg excitation pulses. The ex¬ 
perimental results for the optical density are shown in 
Fig. |6l In order to compare these results to the model 
we calculate values for the decay rates matching the cho¬ 
sen combination of states in our experiment (see Fig. [Ij. 
We calculate the spontaneous emission rates using the 
wavefunctions calculated in [I^ . r 32 is approximated by 
summing the spontaneous decay rates from |35Si/2) to 
all nPj states, which results in 

r 32 = ^ rsp_ 35 S-).nP = 2tt X 3.9 kHz. (12) 

n>5 
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FIG. 7. (Color online) Dynamics of the population of the four 
states (pii=5Si/ 2, p22=5P3/2, P33=35Si/2, p44=Reservoir) 
retrieved from the fit to the dataset with high excitation (see 
Fig. [ 3 . The colors match those of the states in Fig. [3] The 
population of the Rydberg state (35Si/2) starts to decay im¬ 
mediately after excitation at t = 0. 


Rabi frequencies are consistent with estimates based on 
beam power and geometry. The noise 7 c is mainly caused 
by the acoustic noise on the Fabry-Perot cavities. As the 
fit is only sensitive to 7 p -I- 7 c and not to the single values, 
we choose 7 p = 2 7r x 20 kHz. 

An analysis of the datasets with the same conditions, 
but with Rydberg excitation, allows us to fit the frac¬ 
tion of atoms excited to a Rydberg state (0.494(8) and 
0.284(10)) within interval A, where only the red laser is 
on. The superradiance parameter Psnp = 7.9(8) x 10^ is 
fitted to parts B and C and the dipole-dipole interaction 
parameter 73 ,dd is adjusted in part B. The former scales 
with the absolute atom number while the latter scales 
with the atomic density. From the model one can now 
derive the time-resolved populations of the participating 
states as can be seen in Fig. [T] The accuracy of these 
populations is within ±0.01 compared with values calcu¬ 
lated using a variation method. 


VI. CONCLUSION 


The main contribution comes from I 5 P 3 / 2 ) (Fgp, 353 ^. 5 ? = 
27r X 1.2 kHz) and other low-lying, fast-decaying states. 
The spontaneous decay rate r 34 ^sp is given by the rate 
rsp,35S-)'34P = 27r X 16.8 Hz. Here we take only the 
strongest superradiant transition into account. The tran¬ 
sition rate F 34 _bb is approximated by a sum over all tran¬ 
sition rates induced by blackbody radiation from |35Si/2) 
to all nPj states 

r 34 ,bb = ^ Tbb.sss-s-nP = 27r X 2.7kHz, (13) 

n>5 

in which the rates Fbb.sss-s-nP are calculated as in 
The main contribution to r 34 _bb comes from neighboring 
Rydberg states. For the transitions induced by black- 
body radiation a temperature of 300 K is assumed. The 
preceding approximations for the decay rates ensure that 
the total decay out of the state | 35 Si/ 2 ) is modeled cor¬ 
rectly. Similar to the calculation of F 32 we obtain 

r4i = ^ rsp_ 34 P_,.„s = 27r X 0.8 kHz. (14) 

n>5 

The Rabi frequencies of the two lasers Hp = 
27r X 0.83 MHz and He = 27r x 2.10 MHz and the coher¬ 
ence decay 7 c = 2 7r x 112 kHz are htted to the dataset 
without Rydberg excitation (see curve T’ in Fig. ED. The 


We have demonstrated the all-optical detection of Ryd¬ 
berg population in a dilute gas, which is an alternative to 
the methods based on field ionization. Our results show 
that Rydberg population fractions can be measured with 
an accuracy of 0.01. By comparing the dynamics of the 
measured optical densities to our numerical simulations 
we have quantified the decoherence effects occurring in 
the system, namely blackbody radiation induced transi¬ 
tions, superradiant decay and inhomogeneous broadening 
due to dipole-dipole interactions. 

From our simulations we conclude that the detection 
scheme can also be used to obtain information on the co¬ 
herence between the ground state and the Rydberg state. 
The numerical results predict that future studies using a 
coherent excitation method and experimental parameters 
similar to our experiment will be able to detect the initial 
and time-dependent coherence. 
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